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Abstract. This paper is devoted to analyze the null controllability of a 
, nonlinear reaction-diffusion system approximating a parabolic-elliptic system 

modeling electrical activity in the heart. The uniform, with respect to the 
(N ■ degenerating parameter, null controllability of the approximating system by 

means of a single control is shown. The proof is based on the combination of 
Carlemans estimates and weighted energy inequalities. 
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1. Introduction 



Let f2 C R N (N > 1) be a bounded connected open set whose boundary dVl is 
regular enough. Let T > and let uj and O be (small) nonempty subsets of f2, which 
will usually be referred as control domains. We will use the notation Q — f2 x (0, T) 
and E = 50x (0,T). 

In this paper we study the controllability and observability properties for a 
reaction-diffusion system which degenerates into a parabolic-elliptic system de- 
f-C ' scribing the cardiac electric activity in f2 (fi C R 3 being the natural domain of 

the heart). 

\^ ' To state the model, we set = Ui(t,x) and u e = u e (t,x) to represent, respec- 

tively, the spatial cellular and location x £ fl of the intracellular and extracellular 
electric potentials, whose difference v = v(t,x) — Ui — u e is the transmembrane 
potential. The anisotropic properties of the two media are modeled by intracellular 
and extracellular conductivity tensors Mj(a;) and M e (ir). The surface capacitance 
of the membrane is represented by the constant c m > 0. The transmembrane ionic 
current is represented by a nonlinear function h(v) (the most interesting case being 
when the nonlinearity is cubic polynomial) depending on the value of the potential 
v. 

The equations governing the cardiac electric activity are given by the coupled 
reaction-diffusion system: 

Cmdtv — div(Mi(a;)V'Ui) + h(v) = fl w in Q, 
c m d t v + div (M e (x)\7u e ) + h(v) = gl a in Q, 

where / and g are stimulation currents applied, respectively, to cu and O. 
System is known as the bidomain model. 
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We complete the bidomain model with Dirichlet boundary conditions for the 
intra- and extracellular electric potentials 

(1.2) Ui — and u e — on S 
and with initial data for the transmembrane potential 

(1.3) v(0,x)=v o (x), ie!l. 

It is important to point out that realistic models describing electrical activities 
include a system of ODE's for computing the ionic current as a function of the 
transmembrane potential and a serie of additional "gating variables" aiming to 
model the ionic transfer across the cell membrane (see [TTT [T71 12"31 l2~j] V 

Assume fl u = glo- If M.; = ^M e for some constant //el, then system (jl,l|) 
is equivalent to the parabolic-elliptic system: 



{c m d t v - ^div (M e (x)Vv) + h(v) = fl u in Q, 

-div (M(a;)V« e ) = div (Mi(s)V«) in Q, 

v = 0, u e = on S, 

v(0) = vo in il. 

where M — Mi + M e . System (| 1 -4[) is known as the monodomain model. 

We approximate the monodomain model by the family of parabolic systems: 

'c m d t v s - ^div (M e {x)\7vt) + h(v E ) = f s l u in Q, 

ed t u E e - div (M(x)Vit|) = div (M i (a^)V« e ) in Q, 

V = 0, < = on E, 

w £ (0) = u Q , u|(0) = u e ,o in 



Since w = iti — u e in the bidomain model, it is natural decompose the initial 

condition vq as vq = Ui t o — u e ,o- 

The aim of this paper is to give an answer to the following question: 

//, for each e > 0, there exists a control f e that drives the solution (y e ,u e e ) of 

hl.5\) to zero at time t — T , i.e. 

V e (T)=<(T)=0, 

is it true that when e — > the control sequence f e converges to a function f , that 
drives the solution (v, u e ) of Jj.^[ ) to zero at time t = T? 

Since the bidomain model is a system of two coupled parabolic equations and 
the monodomain model is a system of parabolic-elliptic type, these two systems 
have, at least a priori, different control properties. Therefore, it is natural to ask if 
the controllability of the monodomain model can be seen as a limit process of the 
controllability of a family of parabolic systems. 

It is worth to mention that systems of parabolic equations which degenerates into 
parabolic-elliptic ones arise in many areas, such as biology and models describing 
gravitational interaction of particles, see fH 151 12"0]. 
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In order to answer the previous question, we consider the following linearized 
version of (11.51): 



(1.6) 



in Q, 
in Q, 
on E, 
in f2, 



'c m d t v' - -j^div (M e (x)Vv*) + a(t, x)v* = fn u 
ed t u% - div (M(x)Vu e e ) = div (Mi(x)W) 
v e = 0, u e e =0 
K v e (0) = v , u e e (0) = ti e ,o 

where a is a bounded function. 

Our objective then will be drive both v e and u e , solution of (|1.6[) , to zero at time 
T by means of a control / c in such a way that the sequence of controls f e remains 
bounded when e — > 0. Accordingly, we consider the corresponding adjoint system: 

-c m d t (p £ - -j^div (M e (x)Vtp e 
-edt<p e e - div (M(x)Vip%) = 



a(t, x)cp e = div (Mi(x)Vipl 



(1.7) 



<p e = 0, ift = 
^ e {T) = (f T , tp £ e {T) = Lp e , T 



m 

in Q, 
on E, 
in J7. 



It is very easy to prove that our task turns out to be equivalent to the following 
observability inequality: 



(1.8) e||^(0)||£ a(n) + Wm\h (n) ^ C jj Q \^ e \ 2 dxdt, := u x (0, T), 

for all (ipT,<p e ,T) S i 2 (f2) 2 and a constant C = C(e, O, w, T) that remains bounded 
when e — > 0. 

Let us now mention some works that have been devoted to the theoretical and 
numerical study of the bidomain model (|l.lj) . In [6], it is proved the existence 
of weak solutions to (|1.1[) using the theory of evolution variational inequalities in 
Hilbert spaces. Applying the same approach, Sanfelici proves the convergence 
of Galerkin approximations for this model. Recently, Bcndahmane and Karlsen 
PQ proved the existence and uniqueness for a nonlinear version of the bidomain 
equations (|1.1[) by using a uniformly parabolic regularization of the system and the 
Faedo-Galerkin method. 

Regarding finite volume (FV) schemes for cardiac problems, Bendahmane and 
Karlsen [2] analyse a FV method for the bidomain model, supplying various ex- 
istence, uniqueness and convergence results. Bendahmane, Burger and Ruiz [3] 
analyse a parabolic-elliptic system with Neumann boundary conditions, adapting 
the approach in [2]; they also provide numerical experiments. 

Let us now recall some results on the controllability for systems of parabolic 
equations. In [15] the controllability of a quite general two coupled linear parabolic 
system is studied and null controllability is obtained by means of Carleman inqual- 
ities. In [19] , using a different strategy, the controllability of a reaction-diffusion 
system of a simple two coupled parabolic equations is analyzed, the authors prove 
the null controllability for the linear system and the local null controllability of 
the nonlinear system. Another relevant work concerning to the controllability of 
coupled systems is [10] , in which the authors analyze the null controllability of a 
cascade system of m coupled parabolic equations and the authors are able to ob- 
tain null controllability for the cascade system whenever they have a good coupling 
structure. But, unlike the present work, the aforementioned works are devoted to 
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systems that do not degenerate. Actually, if one follow their proofs, it can be seen 
that the constant C appearing in the observability inequality (|1.8|) is of order of e _1 , 
which degenerates when e — > 0. Therefore, a careful analysis is required in order 
to guarantee uniform controllability with respect to the degenerating parameter e. 
This will be done combining precise, with respect to e, Carleman estimates and 
weighted energy inequalities. This kind of analysis has been used several times in 
the case of parabolic equations degenerating into hyperbolic ones (see |16[ l7[ H3]) 
and hyperbolic equations degenerating into parabolic ones (see [HJ[32]) but, as far 
as we know, this is the first time that controllability of parabolic systems degener- 
ating into parabolic-elliptic systems is studied. 

Concerning to the controllability of the bidomain model, since we have couplings 
given by time derivatives of the electrical potentials on both equations, it is very 
difficult to analyze whether the bidomain model is null controllable or not, even 
with controls on both equations. As far as we know, this problem is still open. 
Regarding the monodomain model, since the solution to the parabolic equation 
enters as a source term in the elliptic equation, the following Theorem holds: 

Theorem 1.1. Given vo in L 2 {VL) and 

(1.9) q N G (2, co) if N = 1,2, < q N < 2^±| if N > 3 . 
We have: 

• Ifh is C^R), global lipschitz and satisfies h(0) — 0. Then there exists a control 
f G L 2 (lu x (0,T)) such that the solution (v,u e ) of |_?.^[ ) satisfies 

v(T) = u e (T) = 0. 

• If h is of class C 1 satisfying 

(1.10) mo) = o, h M- h ^ >- C , vt^t*, 

Vi — v 2 

h(v) , h(v) 

(1.11) < liminf — — < limsup — ^— < co, 

\v\-H>o V A |„|->-oa V 6 

andv G H&(n)r\W 2il ~™ ),qN (n), with || VoIIl 00 < 7; for sufficient small 7. There 
exists a control f G L qN (u x (0, T)) such that the solution (v,u e ) of Ji.^l ), with 
(v,u e ) G Wq£(Q) 2 , satisfies 

v(T) = u e (T) = 0. 

Proof of Theorem 11.11 can be obtained from a more general Theorem in [5] . 
For a general discussion about the controllability of parabolic systems see the 
survey paper [18] . 

In this paper we prove null controllability to (|1.6I) for each s > by establishing 
an observability estimate like (|1.8[) for its dual system (|1.7I) . Moreover, the estimate 
on the control we obtain are uniform with respect to e, i.e., the constant C appearing 
in (|1.8I) does not depend on e. In addition, we study the controllability of the 
nonlinear system (|1.5[) obtaining, under some assumptions on the nonlinearity and 
the initial data, uniform null controllability. 

This paper is organized as follows: In Section [21 we state the main results. 
Section [3] is devoted to prove a Carleman inequality for the adjoint system (|1.7[) . 
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In Section |4] we show the uniform null controllability for (|1.6|) . Section [5] deals with 
the uniform null controllability for the nonlinear system (|1.5p . 

2. Main results 

Throughout this paper we will assume that the matrices Mj, j — i,e are C°°, 
bounded, symmetric and positive semidefinite. 
We have the following existence Theorem. 

Theorem 2.1. Under conditions (jl.101) and If (vo,u _ e ) £ L 2 (fl) 2 and 

f £ L 2 (Q), then system \1.5}) have a unique weak solution (v £ ,u £ ) and (v £ ,u £ e ) £ 
L 2 (0,T;H^(n))xL 2 (0,T;H^{n)) such thatd t v e anded t u £ e belong to L 2 (0,T, H 1 ^))- 
L 4/3 (Q) and L 2 (0, T, iJ _1 (f2)). 

The proof of Theorem I2.1l can be done exactly as in [I], so we omit it. 

Our first main result is a uniform Carleman estimate for the adjoint system (| 1 . T|) . 

Theorem 2.2. There exist positive constants C = C(Q,luo), Xq and Sq so that for 
every {pT,^p e ,T £ L 2 (ty, a G L°°{Q), the solution (ip £ ,(p £ e ) of \1. 7| ) satisfies 



e 3sa \p £ \ 2 dxdt + s 3 X 4 // 3 e 3sQ |^ £ | 2 dxd< 

(2.1) < C e 6A||^M s 7 A 4 /" 8 e 2,a|^|2^ 



w/iere // = div (M(x)V^(x, i)), /or every s > (T + (1 + ||a|| 2/ i + ||a||^i)T 2 + 
T )sq, X > Xq and appropriate weight functions <f> and a defined in \S.3\) and \3.$ , 
respectively. 

Proof of Theorem 12.21 is given in Section [3] . 

Our second main result gives the null controllability of (|1.6[) . 

Theorem 2.3. Given vq and w e .o in L 2 (Q). For each e > 0, there exists a control 
f c £ L 2 (u> x (0,T)) so that the associated solution to U.6\) is driven to zero at time 
T . That is, the associated solution satisfies 

v £ (T) = 0, <(T) = 0. 

Moreover, the control f £ satisfies 

(2-2) \\.f e U\mQ } <C{\\vo\\ L2m +s\\u efi \\ L2(n) ). 
Theorem 12.31 is proved in Section 01 

The third main result of this paper is concerned with the uniform null control- 
lability of the nonlinear parabolic system (jl.5l) . 



Theorem 2.4. Given vq and u ei o in L 2 (Q) and let qw satisfying U.9\) . We have: 
• If h is C 1 (M), global lipschitz and satisfies h(0) = 0, there exists a control f £ £ 
L 2 (uj x (0,T)) such that the solution (v £ ,ul) of U.5\) satisfies 

v e (T)=ut(T)=Q. 

Besides, the control f £ has the estimate 

(2-3) \\rU\L HQ) <C(\\v \\ La(n) +e\\u e , \\ L3(sl) ). 
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•Let h be a C 1 function satisfying ALIO)) and 11.11)) and the initial data (vq^u^q) £ 

(Hq(Q) n W 2< - 1 ~«n' I ' 9n (ft)) 2 : with ||(vo,« e ,o)||L°° < 1, for sufficient small 7 does 
not depending on e. There exists a control f £ £ L qN (lo x (0,T)) such that the 
solution {v e ,U%) of 111. 5)) , with {v e ,U%) £ {W 2 ^(Q)) 2 , satisfies 

v s (T) = ul(T)=0. 

Moreover, the control f e has the estimate 

(2-4) \\f e U\lm iQ) <C(|ko|li= ( n)+e||"e,o||" 2(n) ). 

We prove Theorem 12.41 in Section [5l 

3. A Carleman type inequality 



This section is devoted to prove Theorem 12.21 

To simplify the notation, we drop the index e and, since the only constant which 
matters to us is e, we will suppose that all other constants are equal to one. The 
adjoint system (jl.7p then writes as 



(3.1) 



-dt<p — div {M e {x)W<p) + a(x, t)ip = div (Mj(a;)V^e) m Qi 

—edt^Pe ~ div (M(x)\7ip e ) = in Q, 

<f = 0, ip e = on S, 

% <p(T) = tpr, ¥>e(T) = <p e ,T in fi. 



(3.2) 



Suppose that </?t and (/5 e> T are regular enough. Taking p (x, t) — div {Mi{x)Wip e {x, t)), 
we see that the pair (<p, p) satisfies 

—dt<p — div (M e (x)V(p) + a(x, t)(p = p in Q, 

-ed t p - div (M {x)Vp) = in Q, 

(p = 0, p = on S, 

.f( T ) = ¥>T, p{T) = Pt in fi. 

Before start proving the Carleman inequality (|2.ip . let us define several weight 
functions which will be usefull in the sequel. 

Lemma 3.1. Let wo be an arbitrary nonempty open set such that ZJo C w C fi. 
Then there exists a function ip £ C 2 (il) such that 



tj)(x) > 0,Vx £ Q, ijj = 0ondn, \Vip{x)\ > Vx e Sl\u . 
Proof. The proof of this lemma is given in [5]. □ 
From Lemma 13.11 we introduce the weight functions 

e A(V(x)+m|h/>||) e Am|M 

(3.3) <f>(x,t) = — r — ; <j)*(t) = mm(f)(x,t) 



t(T-t) ' r w xe n y ' t(T-ty 
(3.4) 

e A«>(aO+m|h/>||) - e 2Am||V|| gA(m+l)||i/>|| _ e 2Am|h/>|| 



a{x,t) = — — ; a*(t)=maxa(j;,t) 

t(l - t) xe n t(l - t) 

for a parameter A > and a constant m > 1. Here 

HV'Oe)!! = max (V'(^)l ■ 
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Remark 1. From the definition of a and a* we have that 3a* < 2a (for X large 
enough!). Moreover 



4>*{i) < (f){x,t) < e X ^(f>*(x,t) 



\d t a*\<e 2X ^T<f> 2 . 

Now we prove Theorem [ 



Proof of Theorem 
steps:: 



For an easier comprehension, we divide the proof in several 



• Step 1 First estimate for the parabolic system 

In this step we obtain a first Carleman estimate for the adjoint system. 

We consider a set u>\ such that uiq CCwi CCw and apply Carleman inequalities 
(|6.2p . with e = 1, and (|6.15[) to (p and p, respectively, to obtain 



S - 1 (l)- 1 e 2sa \ip t \ 2 dxdt + s~ 



-l e 2sa 



N 



dxdt 



(3.5) +s 3 A 4 fj ^ e 2s a \ v f dxdt + s\ 2 JJ <j>e 2sa \Vip\ 2 dxdt 

: C[ II e 2sa (\p\ 2 + \^\ 2 )dxdt + s 3 X 4 [[ <ft 3 e 2sa \p\ 2 dxdt I , 



and 



(3.6) 



s~ L e 2sa \d t p\ dxdt + e- 



s^e 2 ™ 



N 



dxdt 



+ s 3 \ 4 e- 2 II ^e 2sa \p\ 2 dxdt + s\ 2 e- 2 



2 e 2sa \Vp\ 2 dxdt 



<Ce x ^h 3 X 4 e- 2 // ^e 2sa \p\ 2 dxdt. 



s 
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Next, we add (|3.5p and (|3.6|) and absorb the lower order terms in the right-hand 
side, we get this way 



(3.7) 



L e 2sQ \(p t \ dxdt- 



-l e 2s a 



N 



^ ^ iXj 



dxdt 



+ six4 JJ <f> 3 e 2sa \(p\ 2 dxdt + s 2 \ 2 JJ (j)e 2sa \\7ip\ 2 dxdt 



2sa 



dtp\ dxdt 



2sa 



Q 

N 



dxdt 



+ s 4 A 4 // (f> 4 e 2sa \p\ 2 dxdt + s 2 \ 2 II <j> 2 e 2sa \Vp\ 2 dxdt 



<C[ e x ^s i 



X i JJ 4 e 2sQ |p| 2 dx^ + s 4 A 4 J J <j?e 2sa \^\ 2 dxdt^j. 



for s > (T + (1 + ||a|| 2/ i + \\a\\^l)T 2 + T 4 )s . 

At this point a remark has to be done. If we were trying to control (|1.6[) with 
controls on both equations, inequality (|3.7I) would be sufficient. 



• Step 2 Estimation of the local integral of p. 

In this step we estimate the local integral involving p in the right-hand side 
of (|3.7[) . It will be done using equation (|3.2[h . Indeed, we consider a function £ 
satisfying 



£ G Cg°(w), < £ < 1, £(a;) = lViewi 



and then 

Ce A|Hl s 4 A 4 



~ 2sQ </> 4 H 2 dxdi 



<Ce x ^s A \ i jj e 2sa cj) A \p\ 2 ^dxdt 

= Ce A|Wl s 4 A 4 J J e 2sa c/> 4 p(-<pt-<liv(M e V<p) + acp)tdxdt 
:=E + F + G, 

In the sequel we estimate each parcel in the expression above. 



E = Ce^WjJ sdtae^^p^dxdt 

+ ( 7 e A||Vll s 4 A 4 JJ e 2sa f(j>tp^dxdt + Ce x ^s A X 4 / / , x - < > ^ //V i «'<-<" 
:= i?i + i?2 + £"3. 
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(3.8) 



E 1+ E 2 < 4:S 4 A 4 / / e 2sa ^\p\ 2 dxdt + Ce 2A ^H S 8 A 4 / / e 2sa ^\ 2 dxdt 



10 



and 



E ^yJJ z 2sa \dt P \ 2 dxdt + Ce 2X ^e- 2 s 8 X 8 J f e 2sa cf\ V \ 2 dxdt. 



Next, 



N 

e -A||Vll s -4 A -4 F= V- // s d Xi ae 2sa ((> 4 P {Mi j d Xj ip)tdxdt 



E // 

+ J2 ff e 2sa 4> S d x ^ P {M^d Xj ^dxdt 
+ J2 [f ^ sa ^d x%P {M^d X]V )idxdt 



V // e 2sa cb 4 p(M^d x ^)d x ^dxdt 



and we can show that 



F < ^ S 4 A 4 / 1 e 2s > 4 M 2 ^ + i S 2 A 2 / / e 2s > 2 |Vp| 2 ^ ( /. 



f "z , ■ 

Finally, 



+ Ce 2X ^s 8 \ s JJ e 2sa ^\^\ 2 dxdt+^JJ e 2sa \dl Xj p\ 2 dxdt. 



G<4ts 4 A 4 // e^^lpl 2 ^ 
10 JJQ W 



■ \'\a\U jj " ' - : 
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Putting E, F and G together in Q3.7p . we get 
(3.9) 




Using (|3.9[) we could prove that, for every e > 0, system (11.6[) is null controllable, 
but the sequence of control obtained this way will not be bounded when e — > 0. 
Therefore, we need to go further and improve ()3.9|) . This will be done in the next 
step. 

• Step 3. An energy Inequality. 

The reason why we do not obtain a bounded sequence of controls out of step 2 
is because of the term e~ 2 in the right-hand side of (I3.9|) . In this step we prove 
a weighted energy inequality for equation (|3.2[) 9. which will be used to, somehow, 
compensate this e~ 2 term. 

Let us introduce the function 

V = e^ Q *p, 

which solves the system 

{edty - div (M(x)Vy) = e^sd t a*e? sa ' p in Q, 

y = on E, 

2/(0) = y(T) = in Q. 

We multiply (|3 . 10[) by y and integrate on ft, we get 

~ j t \\y(t)\\h m + C\\Vy(t)\\h m < 4 J n sdta * (t)e^ {t) P (t)y(t)dx. 

Integrating this last inequality form to T and using Poincare's and Young's in- 
equalities, it is not difficult to see that 

(3.11) J J e 3sa *\p\ 2 dxdt<Cs 2 e ix ^ J J s 3 ^e 2sa \p\ 2 dxdt. 
From ([3J]) and (|3~TT|) we obtain 

(3.12) J J e 3sa *\p\ 2 dxdt<Ce GX ^s 7 \ A JJ <P*e 2sa \ V \ 2 dxdt. 

This estimate gives a global estimate of p in terms of a local integral of </>, with 
a bounded constant. 
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• Step 4- Last estimates and conclusion. 

In order to finish the prove of Theorem \2.2\ we combine inequality (|3.12p and an 
another Carleman inequality to equation (|3.2p i . Indeed, we have 



— 1 i — l 3sct 

s m e 



N 



tp t \ 2 dxdt + s" 1 / / ^ 1 e 3sQ | d^ X] (f\ 2 dxdt 



( 3 - 13 ) + S 3 \*JJ < t )3e3Sa ^\ 2 dxdt + s\ 2 JJ <j)e 3sa \Vip\ 2 dxdt 

~ C {.II e3Sa \ p \ 2dxdt + s3X4 JJ < t )Ze3sa \ i P\ 2dxd ^ 

where <p is, together with p, solution of (|3.2j) . 

Here, we just changed the weight e 2sa by e 3sQ . The proof of (|3.13|) is the same 
as the one given by Theorem 16. 11 just taking a slightly different change of variable 
in (IO|) . 

Next, since e 3sQ < e 3sa , we have 

„3sa\ „|2 j„ j + II „3sa*|„|2 



e^H^dt < / / e isa \p\ 2 dxdt 



and by (j3~T2]) . 



e 3sQ |p| 2 da;di < Ce 6A|wl s 7 A 4 / / 8 e 2sQ |^| 2 dxdi. 

From (|3T2|) and f3l3]> . it follows that 
(3.14) 

/ / e 3sa \p\ 2 dxdt + s 3 \ 4 f f <f> 3 e 3sa \cp\ 2 dxdt < Ce 6X ^s 7 X 4 f f <fe 2sa \^\ 2 dxdt, 

JJq JJq JJq^ 

which is exactly (|2.1[) . 

By density, we can show that (|3.14j) remains true when we consider initial data in 
L 2 (fl). Therefore, the Carleman inequality (12.11) holds for all initial data in L 2 (fl). 

This finishes the proof of Theorem 12.21 

□ 

4. Null controllability for the linearized system 

This section is intended to prove the null controllability of linearized equation 
(|1.6[) . It will be done by showing observability inequality (|1.8[) for the adjoint system 
(|1.7p and solving a minimization problem. The arguments used here are classical 
in control theory for linear PDE's, so that we just give a sketch of the proof. 



• Sketch of the proof of Theorem \2.3i By standard energy inequality for system 
|), we can prove, using (|3.14[) , that 



(4-1) e||p(0)||l 2(f2) + ||^(0)||| 2W <C ie c ^ T JJ Q \^\ 2 dxdt 

for some constants Ci, C2 > 0. 

Since p(x,t) — div (M(x)Vip e (x, t)) and ip e = on <9f2, we have 

IM*)||jr»(n) <C\\p(t)\\ L2m , 
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for all te[0,T]. Therefore, it follows from |@~TJ that 




which is the observability inequality (|1.8[) . 

From (|4.2[) and the density of smooth solutions in the space of solutions of (|3.1[) 
with initial data in L 2 (f2), we see that the above observability inequality is satisfied 
by all solutions of (jl.7j) with initial data in L 2 (fl). 

Now, in order to obtain the null controllability for linear system (|1.6j) . we need 
to solve the following minimization problem: 

Given ipT and <p e ,T in L 2 (il), 

Minimize J'g (ifT , <Pe,T)i with 

J S (<PT,<Pe,T) = \~ [ [ \f £ \ 2 dxdt + £(Ue,Q,<p S e (0)) 

(4.3) [ 2 Jo Jcu 

+ (V ,<p e (0)) + 6(\\<pT\\ LHQ) + |ke,T|| i2(n) )|, 

where (tp, tp e ) is the solution of the adjoint problem (|1.7|) with initital data (yr, V?e,T)- 
It is an easy matter to check that is strictly convex and continuous. So, in 
order to guarantee the existence of a minimizer, the only thing remaining to prove 
is the coercivity of J$ . 

Using the observability inequality (|1.8j) for the adjoint system (jl.7j) . the coer- 
civity of Js is straightfoward. Therefore, for each 5 > 0, there exists a unique 
minimizer (ip s eT ,ip^) of Js. Let us denote by ip £ ' S the corresponding solution to 
(|1.7p associated to this minimizer. Taking f £ - s = ip s ' 5 l u as a control for (|1.6|) . the 
duality between (|1.6j) and (|1.7j) gives 

(4-4) \\v s > s (T)\\ L i m + \\uf(T)\\ L2{n) <S, 

where (v e ' S ,u £ e ' S ) is the solution associated to the control f e ' S . It also gives 

(4-5) \\f s ' 5l u\\L2(Q) < C(||w |U=(n) +s\\u ei0 \\ L 2 {u) ). 

From (|4.4p and (|4.5I) . we get a control f e (the weak limit of a subsequence of f e ' 5 lui 
in i 2 (w x (0,T))) that drives the solution of (|1.6[) to zero at time T. From (|4.5j) 
we have the following estimate on the control / e , 

(4-6) \\f e U\LHQ) < C(H\\ L2m +s\\ul fi \\ LHQ) ). 

This finishes the proof of Theorem 12.31 □ 



5. The nonlinear system 

In this section we prove Theorem 12.41 The proof is done applying fixed point 
arguments. 
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Proof of Theorem \2.J\ ( case 1 ): We consider the following linearization of system 
LB- 

'c m d t v E - ^div (M e (x)Vv E ) +g{z)v s = f E l u in Q, 

ed t u e e - div (M(x)Vul) = div (M z (x)\7v £ ) in Q, 

v e = 0, u% = on E, 

k u e (0) = v , u|(0) = it e ,o in Q, 



(5.1) 



where 

(5.2) g( S ) 



^ if \s\ > 0, 

s 

ft'(s) if s = 0. 



Follows from Theorem 12 .31 that for each «o, «e,o 6 L 2 (f2) and z 6 L 2 (Q), there exist 
a control function f e Xui £ L 2 (Q) such that the solution of (|5.ip stisfies 

V £ (T)=^(T) = 0. 

As we said before, the idea is to use a fixed point argument. For that, we will use the 
following generalized version of Kakutani's fixed point theorem, due to Glicksberg 

Theorem 5.1. Let B be a non-empty convex, compact subset of a locally convex 
topological vector space X . If A : B — > X is a set-valued mapping convex, compact 
and with closed graph. Then the set of fixed points of A is non-empty and compact. 

In order to apply Glicksberg's Theorem, we define a mapping A : B — > X as 
follows 

A(z) = {v; (v,u e ) is a solution of (|5.ip . such that v(T) = u e (T) = 0, 

for a control fxu satisfying (|2.2[) }. 

where X — L 2 (Q) and B is the ball 

B = {z e L 2 (0,T,H^(n),d t z S L^O.^ff-^Q)); 

W z \\ 2 L 2 (o,T;H},(n)) + H^lli^o.Tjff-^n)) - ^0- 

It is easy to see that A is well defined and that B is convex and compact in L 2 (Q). 
Now, we prove that A is convex, compact and has closed graphic. It will be done 
into the next steps. 

• A(B) C B. 

Let z G B and v G A(z). Since v satisfies (|5.ip i . the following inequality holds 

( 5 - 3 ) \\ V \\i 2 (Q,T-Hi(n)) + ll^t U lli a (Q,T;H- 1 (n)) - K ^ 

In this way, if z £ B then A(z) C B, if we take M = K\. 

• A(z) is closed in L 2 (Q). 

Let z G B fixed, and v n G A(z), such that v„ — > v. Let's prove that v G A(z). 
In fact, by definition we have that v n is, together with a function it e! „ and a 
control /„, the solution of (I5TT1) . with ||/ n lo)|U=(Q) < C(|KIU2(n) +£||w e ,o|U 2 (n))- 
Therefore we can extract a subsequence of /„, denoted by the same index, such 
that 

fn^-u, -> /1« weakly in L 2 (Q). 
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Since f n is bounded, we can argue as in the previous section in order to obtain the 
inequality 

(5-4) H U nlL 2 (0,T;Hj(n)) + \9tV n | \ L 2 ( 0tT;II -i (O)) ^ @ 

and follows that 

v n v weakly in L 2 (0,T; Hq(SI)), 
v n — > v strongly in L 2 (Q), 
d t v n -> <9 t v weakly in L 2 (0, T; ff" 1 ^)). 

Using the converges above and (|5.1[) 9. we see that there exists a function it e such 
that 

u e ,n -> u e weakly in L 2 (0, T; if* (0)), 
Me, n -> u e strongly in L 2 (Q), 
d t u e<n -> 5 tMe weakly in L 2 (0, T; 

It is immediate that (it e , i>) is a controlled solution of (|5.1I) associated to the control 
/. Hence v G A(z) and A(z) is closed and compact of L 2 (Q). 

• A has closed graph in L 2 (Q) x L 2 (Q). 

We need to prove that if z n — > z, v n — > v strongly in L 2 (Q) and v n G A(z n ), 
then v G A(z). Using previous steps, it is straightforward that v G A(z). 

Therefore, we can apply Glicksberg Theorem to conclude that A has a fixed 
point. This proves Theorem 12.41 in the case which the nonlinearity is a C 1 global 
Lipschitz function. 

Proof of Theorem \2.4\ (case 2): The proof of the local null controllability in the 
case 2 is done exactly as in the equivalent one in |19j . 
We consider the linearization 



(5.5) 



'c m d t v s - ^div (M e (a;)Vw £ ) + a{z)v £ = f" X u in Q, 

ed t u £ e - div (M(x)Vu|) = div (M;(x)W) in Q, 

v £ = 0, u% = on E, 

■y £ (0) = vq, u%(0) = w e ,o in f2, 



with («d,« e , ) G (ff 1 (fi)nlf 2(1_ wi' w (n)) 2 ) z G and 

a(z) = / —(sz)ds. 
Jo dz 

It is not difficult to show the null controllability of (| 1 . 5|) with a control in L 2 (cj x 
(0, T)), but these kind of controls are not sufficient to use fixed point arguments in 
order to control the nonlinear system (|1.4p . Our strategy then will be to change a 
bit the functional (|4.3p in order to get controls in L qN (Q) and then apply a fixed 
point argument. In fact, we consider the functional: 

Minimize J'g (ifT , ¥>eT ) , with 

T r 

2sa 




Js(<PT,<PeT) Ho/ I e"Wr dxdt + E(u efi ip e e (0)) 

(5.6) 



(vo,<^ e (0)) + s(\\ip T \\ L 2 {n) + ||v e ,T|| L2(n) 
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where (ip E , tpf,) is the solution of the adjoint system (|1.7|) with initital data (ifT, Ve.r)- 

As before, it can be proved that t|5 ,6|) has an unique minimizer (<p £ ' S , 
Defining f £ - s — e 2sa (f> s (p e ' S and using the fact that <p e ' S is, together with a ip e e s , 
the solution of (11.71) . we see that f £ ' S is a solution of a heat equation with null 
initial data and right-hand side in L 2 (Q). Using the regularizing effect of the heat 
equation and arguing exactly as in |19j we can show the estimate 

(5-7) \\f eJ U\L«N( Q ) <C(\\v f L 

Taking the limit when 5 — > 0, we get a control / e l w G L qN (Q) such that the 
associated solution (v E ,u e e ) to (|5.5j) satisfies 

« E (T) = <(T) = 0. 

The proof is finished by applying Kakutani's fixed point Theorem for system (|5.5|) . 
exactly as done in Theorem 6 in |19j . 

6. Appendix: Some technical results 
Let g G L 2 {Q) and vt G L 2 (J7). In this section we will consider the equation 

N 

-d t v(x,t) - ^ 9 xi(aij(x)d Xi v(t,x)) = g(x,t) in Q, 

,3.1) ^ l -> =1 

v = on E, 

w(T) = w-r in fi, 

under the assumption that the matrix <Zy has the form: 

_ ^% 

where (My)y is an elliptic matrix, i.e., there exists j3 > such that J^fj > 
/3|£| 2 for all £ G R^. 

6.1. A degenerating Carleman estimate. 

Theorem 6.1. There exists Ao > 1 and sq such that, for each, A > Ao and s > 
sq(T + T 2 + T ) i/ie solution v of the equation &6.1\) satisfies 

[[ s- 1 (t)- 1 e 2sa \dtv\ 2 dxdt + s^e- 2 [[ ^- 1 e 2sQ | V dl iX .v\ 2 dxdt 

JJq JJq iJ=1 1 3 

(6-2) +s 3 X 4 e- 2 JJ 4> 3 e 2sa \v\ 2 dxdt + s\ 2 e~ 2 J J 4>e 2sa \Vv\ 2 dxdt 

e2Sa \9\ 2dxdt + s3x4£ ~ 2 ff <t> 3 e 2sa \v\ 2 dxd?J, 
with C > depending on fl, ujq, ip and /3. 
Proof. For s > and A > we consider the change of variable 
(6.3) w(t,w) = e sa v{t lW ), 

which implies 

w(T, x) — u)(0, x) = 0. 
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We have 

(6.4) L 1 w + L 2 w = g s , 
where 

N N 

(6.5) L\w — —dtw + 2sX 4>aijd Xj ip d Xi w + 2sX 2 <f>aijd Xi ip d Xj ipw, 

W AT 

(6.6) L 2 w — - 2J d Xz {a,ijd Xj w) - s 2 X 2 2J <j> 2 aijd Xi ip d Xj ip w + sd t aw, 

»>j=i «ii=i 

and 

JV A 

(6.7) ,g s = e s "g + sA 2 (j>aijd Xi tfj d Xj ipw - sX ^ 4^x t (aijd Xj ip) w. 
From dO), 



(6.8) ||^i^|||2(q) + ||i2-w||| 2(Q) +2(Liw, J L 2 w) L 2 (Q) = ||.9 s ||i2 (Q) . 



The idea is to analyze the terms appearing in {Liw, L 2 w)i,2<q\. First, we write 

A 



(L 1 w,L 2 w) L 2 {Q) = h 
1,3=1 



where Jjj is the inner product in L (Q) of the ith term in the expression of L\w 
and the jift. term of L 2 w and, after a long, but straightforward, calculation, we can 
show that 



2{L lW ,L 2 w) L 2 {Q) >2s 3 X 4 p 2 s- 2 JJ (/> 3 \Vip\ 4 \w\ 2 dxdt + 2sX 2 l3 2 E- 2 JJ (j)\Vip\ 2 \Vw\ 2 dxdt 

-Cs- 2 (t a s 2 X a + Ts 2 X 2 +T 2 s + s 3 X 3 +Ts 2 X^j J J <j> 3 \w\ 2 dxdt 

(6.9) -Ce- 2 (sA + A 2 ) J J cj)\Vw\ 2 dxdt. 

We take A > A and s > s (T + T 2 + T 4 ), it follows from Remark H that 



2{L lW ,L 2 w) L 2 {Q) + 2s 3 X 4 p 2 e- 2 / / 4> 3 \w\ 2 dxdt 



+ 2sX 2 p 2 e- 2 / / (j)\Vw\ 2 dxdt 
(6.10) > 2s 3 X i /3 2 e' 2 J J <f> 3 \w\ 2 dxdt + 2sA 2 /3 2 eT 2 J J <f)\Vw\ 2 dxdt. 



Remark 2. Since f2\o;o is compact and |V"0| > on f2\wo, there exists 8 > such 
that 



/3|V^| > 8 on Q\lu q . 
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Putting ([QUI in (|?T8j) . we get 

\\Liw\\h(Q) + \\L2w\\ 2 L2(Q) +2f3- 2 s 3 \ i 6 i e- 2 jf <jy 3 \w\ 2 dxdt 



3.11) 



+ 2s\ 2 5 2 e- 2 // d>\Vw\ 2 dxdt 



< \ \9s\\h iQ) +2r 2 s 3 X 4 S 4 s- 2 II 3 \w\ 2 dxdt 
+ 2s\ 2 6 2 e- 2 II <b\Vw\ 2 dxdt. 



Now we want to deal with the local integral involving Vw in the right-hand side 
of (|6.11j) . To this end we introduce a cutt-off function £ such that 

£ e C °°H, < £ < 1, t(x) = 1 Va; £ w - 

Utsing the ellipticity condition on , we prove that 

^ II ^\ Vw \ 2dxdt 

< C (^JJ L 2 wcj)t 2 wdxdt + (sT + e-\s 2 \ 2 ) JJ cj) 3 \w\ 2 dxdt 

+ Xe- 1 JJ ^^Vw^^wdxdt). 
By Young's inequality we have 



sX 2 S 2 e~ 2 JJ (f>£ 2 \Vw\ 2 dxdt 

<\JJ \L 2 w\ 2 dxdt + Cr 2 s 3 X i (5 i + S 2 )e~ 2 J j <f) 3 \w\ 2 dxdt. 

Thus, inequality (|6.11|) gives 

\\Liw\\ 2 L2(Q) + \\L 2 w\\ 2 l2{Qt) + /TVaV 2 jj <p 3 \w\ 2 dxdt 
(6.12) +s\ 2 e~ 2 jj (j)\X7w\ 2 dxdt 



Q 



<C[\\e sa g\\ 2 L2(Q) + fi- 2 s 3 \ i e~ 2 // (j) 3 \w\ 2 dxdt 



V 



Now we use the first two terms in left-hand side of (|6.12j) in order to add the 
integrals of |Aw| 2 and \wt\ 2 to the left-hand side of (|6.12|) . This can be made using 
the expressions of L\w and L 2 w. Indeed, from (|6.5|) and (|6.6|) . we have 
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N 

s-^-^dtwfdxdt + e- 2 // s- 1 ^' 1 V \d Xi (Mijd Xj w)\ 2 dxdt 



Q 



(6-13) +. S 3 A 4 £- 2 // <p 3 \w\ 2 dxdt + sX 2 s- 2 // cp\Vw\ 2 dxdt 



< Cy\e aa g\\\ HQ) +s 3 \ i e- 2 JJ <p 3 \w\ 2 dxdtj . 

Using the term in \d Xi (Mijd Xj w)\ 2 in the lef-hand side of (|6. 13|) and elliptic regu- 
larity, it is easy to show that 

s^e- 2 JJ 4>~ 1 \Y1 d l t x,M 2 dxdt < c(\\e sa g\\ 2 L2{Q) + s 3 \ 4 e- 2 JJ (t> 3 \w\ 2 dxdtj . 
Estimate (|6.13[) then gives 



+s 3 X' l e 



s - 1 <j)- 1 \d t w\ 2 dxdt + s^e- 2 { { (f)- 1 \ dl zXj w\ 2 dxdt 
J J (j> 3 \w\ 2 dxdt + s\ 2 e~ 2 JJ <j)\Vw\ 2 dxdt 



3-v4„-2 



' Q 

(6.14) <C[\\e sa g\\ 2 LHQ) +s 3 X 4 e- 2 1 1 ct) 3 \w\ 2 dxdt 



From (|6.14p and the fact that w = e sa v, it is not difficult to finish the proof of 
Theorem O □ 



6.2. A Slightly changed Carleman inequality. 

Theorem 6.2. There exists Xq > 1 and sq such that, for each, X > Xq and s > 
sq(T + T 2 + T ) the solution v of the equation 116.1)) satisfies 

s- 1 e > ' a \dtv\ a dxdt + e- a [[ s^e 2 ™] V d 2 x . x .v\ 2 dxdt 

Q JjQ 4jJ . =1 ' 3 

(6-15) +S ^x i e- 2 JJ ^e 2sa \v\ 2 dxdt + sX 2 e- 2 JJ <t> 2 e 2sa \^ p\ 2 dxdt 

<Ce x ^([[ 0e 2sa \g\ 2 dxdt + s 3 X 4 e- 2 [[ ^e 2sa \v\ 2 dxdt), 



with C > depending on Q, u)q, ip and f3. 

Proof. The starting point is to apply the Carleman inequality given in Theorem 
16.11 for v, i.e., 
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N 

e 2 J J s- 1 q)- 1 e 2sa \v t \ 2 dxdt + J J s-V _1 e 2sQ | 5^ d 2 XzXj v\ 2 dxdt 

( 6 - 16 ) +s3A4 // ^ e2Sa \ v \ 2dxdt + sX2 JJ ^ 2sa \Vv\ 2 dxdt 

<C(JJ e 2sa \g\ 2 dxdt + s 3 X 4 J J (f> 3 e 2sa \v\ 2 dxdt). 
Next, we introduce the function y (x, t) — v(x,t)((f>* (t))^ , which solves the system 



(6.17) 



ed tV - div {M(x)Vy) = -e^Vy + (4>*(t))$g in Q, 

y = on E. 

Applying again the Carleman inequality given by Theorem 16. 11 at this time for y, 
we obtain, for s large enough, that 

s-^^e^ldtyfdxdt + e- 2 [ [ s' 1 cp- 1 e 2sa \ V d 2 x y\ 2 dxdt 
( 6 - 18 ) + ,s 3 AV 2 / / <jy 3 e 2sa \y\ 2 dxdt + sAV 2 [ [ ^e 2sa \Vy\ 2 dxdt 



Q 



<C(j I (f>*e 2sa \g\ 2 dxdt + s 3 A 4 e~ 2 J J 3 e 2sa \y\ 2 dxdt) . 



From the definition of y it is easy to show that 



s -i0-i e 2s > t (0*)2| 2 (fed* < // s- 1 4>- 1 e 2sa \d t y\ 2 dxdt+ // e 2sa cP\y\ 2 dxdt 

Q JJQ JJQ 

and inequality (|6.18[) becomes 
(6.19) 

N 

s^^^e^lvtfdxdt + e- 2 // ^V^V^I ^ d 2 iXj v\ 2 dxdt 

+ s 3 X 4 e- 2 JJ (t> 3 (t)*e 2sa \v\ 2 dxdt + s\ 2 e- 2 JJ 00* e 2sa \Vv\ 2 dxdt 

<C{JJ 0*e 2sa \g\ 2 dxdt + s 3 \ 4 e~ 2 JJ cfcp*e 2sa \v\ 2 dxdt). 
From Remark [1] the result follows. □ 
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